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The  Buckllnp^  of  a  Thick  Circular  Plate 
Ucin^  a  Hon-Liiiear  Theory 

Chester  B.  Sensenig 


1 ,   Introduction 

The  object  of  this  paper  is  to  study  the  stability  of  equllibriuir 
of  an  isotropic  circular  cylinder  iJhen  compressed  along  its  curved 
lateral  surface.  The  compression  on  the  curved  part  of  the  boundary 
is  assumed  to  be  such  that  the  generators  of  the  cylinder  remain  ver- 
tical straight  lines,  but  so  that  no  shear  stress  is  developed:  it  is 
as  though  the  cylinder  x^^ere  compressed  by  shrinking  a  very  stiff 
greased  ring  down  on  it.   The  plane  faces  are  assumed  to  be  free  of 
stress. 

The  problem  is  treated  by  m.aking  use  of  a  non-linear  three-dimen- 
sional theory  of  elasticity,  of  the  same  sort  as  treated  in  the  book 
of  Green  and  Zerna.  The  formulation  of  the  theory  used  in  this  report 
follows  that  of  Fritz  John  [1]  in  using  Lagrange  coordinates  and  a 
matrix  notation.  The  essential  point  is  that  no  restrictive  assump- 
tions are  made  concerning  the  thickness  of  the  plate:  the  problem  is 
treated  as  three-dimensional  throughout.   The  displacements  and 
strains  are  not  assiomed  to  be  small.   The  strain  energy  function,  from 
^^rhich  the  stress-strain  relations  come,  could  in  principle  be  chosen 
arbitrarily  without  causing  any  essential  difficulties;  iN^ithin  this 
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report  a  special,  but  for  raany  materials  reasonable,  choice  for  the 
strain  energy  function  is  made. 

When  the  cylinder  is  corapressed  under  the  conditions  stated 
above,  a  very  simple  state  of  strain,  i.e.,  uniform  plane  strain, 
satisfies  all  conditions  of  the  exact  non-linear  theory.  A  small  per- 
turbation with  respect  to  this  simple  state  is  then  made  and  the  ques- 
tion is  raised  as  to  whether  other  possible  equilibrium  states  - 
buckled  states,  in  other  words  -  exist  within  the  framework  of  the 
basic  three-dimensional  theory,  and  vjhat  pressure  is  needed  to  cause 
them.  Such  critical  pressures  are  found,  in  the  present  case  (infi- 
nitely many  of  them  in  fact),  together  with  the  corresponding  modes 
of  bucklingo  This  report  is  restricted,  however,  to  buckled  modes 
with  symmetry  about  the  axis  of  the  cylinder.   For  these  modes  the 
critical  buckling  pressures  and  modes  are  given.   They  are  the  same 
as  those  furnished  by  the  classical  thin  plate  theory  in  the  limit 
case  of  small  thickness  of  the  cylinder. 

It  should  be  pointed  out  specifically  that  these  results  are  dif- 
ferent from  those  given  in  an  earlier  report  by  S.  Lubkin  [2],   Lubkin 
treats  the  buckling  of  a  hollow  cylinder  (as  well  as  the  hollow  sphere 
and  the  rectangular  solid),  but  he  considers  buckling  modes  in  the 
plane  of  the  faces  of  the  cylinder,  i.e.,  cases  in  which  all  quanti- 
ties are  independent  of  the  coordinate  along  the  axis  of  the  cylinder. 

2 •   Brief  Description  of  the  Non-Linear  Theory 
Let  x'  =  (x'  xJ,x')  and  x  =  {x^,x^,x^)   be  the  position  vectors 
to  the  same  material  particle  in  the  strained  and  unstrained  states 
respectively.  The  Lagrange  coordinates,  x.(i  =  1,2,3) >  are  used 
throughout  this  paper  as  independent  variables,  and  the  Eulerian 


3 
coordinates  x!  (1  =  1,2,3)  a^c  alx-;ays  used  as  dependent  variables. 

Let  p..  =3x!//j)r..   (i,j  =  1,2,3),  and  let  p  be  the  matrix  (Pj^J« 
1 J     1    J  ■'■J 

In  the  theory  which  we  are  using  the  raaterial  is  described  by  a  strain 

energy  function  u  i^ihich  is  a  function  of  p« 

Let  q.  .  =9u/5p.  .  (i,j  =  1,2,3),  and  let  q  be  the  matrix  (q.J. 

If  vje  neglect  body  forces,  the  equilibrium  equations  in  Lagrange  form 

are  q.  ..  .  =  0  (i  =  1,2,3)  v:here  (  ),.=?>{    )/^x.,  and  we  will  use  the 
i  J ,  J  J  J 

usual  summation  convention  throughout. 

The  stress  matrix  'u  is  given  by  '"d  =  l/|p!  qp"  where  |p|  is  the 
determinant  of  p  and  -;:-  denotes  the  transpose. 

Let  S  and  S '  denote  the  same  material  surface  in  the  unstrained 
and  strained  states  respectively.  Let  ^  and  5'  be  unit  outer  normal 
vectors  to  S  and  S'  respectively.   Then  t^.^^ldS'  =  q.  .^.dS  where  dS « 

■^  J   J  0       0 

and  df  are  the  elements  of  area  of  S'  and  S  respectively. 

3 •   Strain  Ener,s:y  Function 

In  order  to  apply  the  theory  mentioned  x-jq   specify  the  strain 
energy  function  to  be  used.  VJe  use  the  usual  strain  energy  function 
of  the  linear  theory,  namely  U  -  (A/2)[e]   +  [x[e  ]  'Jhere  A  and  |j,  are 
the  lame  constants,  e  is  the  strain  matrix,  and  [a]  is  the  trace  of 
the  matrix  a.   Since  we  want  U  to  be  a  function  of  p,  we  next  define 
e  as  a  function  of  p  using  a  definition  due  to   K  .  0.  Friedrichs. 
Let  c  be  the  unique  orthogonal  matrix  for  i-Jhich  cp  is  syrrimetric  and 
positive  definite,  and  let  e  =  cp  -  1.  VJe  think  of  the  rotation  c 
as  the  rotational  part  of  the  deformation. 

Next  we  i^rlsh  to  compute  the  matrix  q.  VJe  observe  that  (e  +  1) 
=  p''*p,  e^  =  p'''p  -  2e  -  1,  [e^]  =  [p'''p]  -  2[e]  -  3,  and  U  =  (V2)[eJ^ 
-  2M-[e]  +  M-tp""'?]  -  3m..   Hence  q  =  (A[e]  -  2|i)  (d[e j/3p^j^)  +  2up. 


s, ,    . 


i    t 


Next  we  shox^  that  (a[ e ]/;>p^j^)  =  c'\    ^  ''      From  cc"  =  1  and  cp  =p"c" 
we  obtain  9c''y3P4i,  =  -  c'"'i'Sc/f^-p.^,)c''\   and  hence 

i^.   p]  =  [p^^  ^^]  =  -  [p%^^  ^  c^^J  =^  -  [cp  ^-  en   =  -   [^^   c%p] 
^Pik         ^P-ik  '^Plk  3Pxk         '^'Plk 

=  -  [:^:rr5 —  p] 

^^Ik 

so  that  [3c /ap^^  p]  =  Oc  Thus  ^[e]/o»p^j^  =  3[cp]/9p^^  =  [-)c/«P3_^  p] 

+  [c(ap/Ap^,^)]  =  [c(9p/?5P^j^)]  =  Cj^^. .  We  now  have  q  =  (A[e]  -  2[i)c'"'+2[i-p, 

We  vjill  also  be  vjorking  vjith  linear  combinations  of  derivatives 
with  respect  to  the  p.,  's.  For  this  we  introduce  a  notation  due  to 
Fritz  John  [1],  For  an  arbitrary  3x3  real  scalar  raatrix  c,   let 
I)    (    )    =(',)(  )/-:^-p.,  )^.,  where  (  )  is  either  a  scalar  function  or  matrix 
function  of  the  P^^i^'s. 

We  now  compute  D^q.-,  for  arbitrary  ^  evaluated  for  p  equal  a  dia- 

gonal  raatrix  i-jlth  positive  diagonal  elements.  For  any  p  vie   have 

D^[e]  ^  (^[e]A)p.  .)?.  .  =  c.,^..  =  [c^J,  When  p  is  diagonal  with  posi- 
c,  •'•J-'-J    J-'--'-J 

tive  diagonal  elements,  then  c  =  1  and  hence 

Prom  cc'"  =  1  and  cp  =  p'"'c'"'  ^^^e  obtain  D^c  +  D,_c"  =  0  and  (D^c)p  +  £, 
=  ^"''''  +  p(D^c'"")  for  the  p  under  consideration.  Hence  for  such  p  we 
have  (D  c)p  +  p(D.c)  =  ^"'  -  ?,  and  D^c^^_,  =  Ci^i-^ik/Pii+Pki.* 


1 

The  derivation  given  here  is  due  to  Fritz  John, 


i  .. 


0  (iK:; 


'   1"'  ".''  f  '"   •"' 


,  i-.  . 


a^'.::".    L.!;v 


,  I  ■  .'  '.  V  . 


We  now  have 

D,q.,  =  A5.,  E..+   2li.?.,  +  (Ap  .  .  -  3A  -  2(i)      ^ — ^   - 
4^ik     Ik'^jj    ^-ik    ^    '  33        -^  ^     Pii  +  P^ij   ' 

for  p  diagonal  with  positive  diagonal  elements, 

Belov;  v/e   list   and  niunber   the   results   obtained   thus   far, 

(3.1)  P  =   (Pi-.)   =   Ox!/ax   ) 

(3.2)  e  =  cp   -  1,    c    is   orthogonal,    cp   is   syrrimetric   and   positive 
definite 

(3.3)  U  =   (A/2)[e]2   +  tx[e^] 
(3.1+)        q   =    (A[e]    -  2ia)c'"'  +  2ixp 

(3»5)        q.-  •    .   =  0    (i  =  1,2,3)    fo^   equilibriura  with  zero  body  forces 

(3.6)  Z=    (l/|p|)qp''- 

(3.7)  U^j.CldS'    =  QijgjdS    (i   =   1,2,3) 

(3.8)  D^q,^  =  A5.^^  .  .  .  2,q^  +   (Ap.  .   -  3X  -  2,)   ^i^^-gi 

for  p  diagonal  with  positive  diagonal  elements, 

I4.,   Statement  of  the  Problem  and  Trivial  Solution 
VJe  are  nov;  in  a  position  to  give  a  precise  formulation  of  our 
problem.  As  before  we  let  5'  =  (^j,^^,^')  be  the  unit  outvjard  normal 
vector  to  a  material  surface  in  the  strained  state,  and  we  let 
E,  =   (5-,,^p,^-)  be  the  unit  outward  normal  vector  to  the  same  material 
surface  in  the  unstrained  state.   In  Lagrange  coordinates  we  wish  to 
solve  the  equilibrium  equations  q,.  .  =  0(1  =  1,2,3)  subject  to  the 
follovjing  boundary  conditions: 


U»L) 


X. 


V 


R 


/■■ 


I 


/ 


/ 


21i 
A. 


PIGUP£   Is        Unrestrained   State 


X, 


X. 


PIGUPlE  II;   First  Buckled  Kode 
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(i;,l)   x'  =  (1  -  e)x^  ,  x'  =  (1  -  s)Xp,  and  '"0*0 --'? '•  =  0  on  the  surface 
X,  =  R  cos  Q   and  Xp  =  R  sin  Q  where  e  is  an  arbitrary  con- 
stant such  that  0  <  s  <  1, 

iko^)       T;..  .-C'  =  0  (i  =  1,2,3)  on  the  surfaces  x_  =  +  h,. 

Using  (3»7)  we  see  that  the  boundary  conditions  can  be  expressed 
as  folloi'js: 


(Ij-.B)   xj_  =  (1  -  e)x^,  x^  =  (1  -  e)x2,  and  Qoj^j  =  q^iCos  0  +  q    . 

sin  @  =  0  for  x^  =  R  cos  Q  and  Xp  =  R  sin  Q, 
(I1..I+)    q^.£,^   =  +  q^^  =  0  (i  =  1,2,3)  for  x^  =  +  h. 


X 


First  we  look  for  a  trivial  solution  of  the  form  x'  =  6,x. 

i    11 

(i  =  1,2,3)  where  each  5.  is  a  positive  constant,  ?or  such  a  solu- 
tion p  is  a  constant  diagonal  matrix  with  positive  diagonal  elements. 
Hence  c  =  1  and  by  {3»k)    q  is  a  constant  diagonal  matrix.   Hence  the 
equilibriiun  equations  are  satisfied,  q,^  cos  Q  +  q-p  sin  Q  =  0  for 

,  =  R  cos  Q  and  x„  =  R  sin  ^,  and  q-, -,  =   q^^  =  0  for  x-,  =  +  h.  From 
the  other  boundary  conditions  we  obtain  5,  =  5p  =  1  -  e  and  0  =  q_^ 
=  A(5^  +  5p  +  5  )  -  3A  -  Z\\.   +  2iJ,5-,   Therefore  5^  =  1  +  (2Ae/X+2n.), 
and  our  trivial  solution  is 


X.  =  (1  -  e)x^  ,   x^  =  (1  -  £)X2   ,   xt  =  (1  +  |^)X3 


5^   Buckled  Solutions 

VJe  noi-J  assume    the    existence   of   a   solution  in  vjhich  x'    (i   =   1,2,3) 

i 

and  e  depend  analytically  on  a  parameter  5  in  the  following  way: 

o       •  •  ' 

x'  =x'  +x'5  +  e..   and   e  =  e  +  e5  +  ... 
1    i    i  o 
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o 

where  x|  (i  =  1,2,3)  is  tho  trivial  colutlon  corresponding  to  e  =  e 

and  where  xj  (1  =  1,2,3)  is  non-trivial  for  5  >  0.   Then  xj  (  i  = 

1,2,3)  is  a  trivial  solution  for  5=0  and  a  buckled  solution  for 

5  >  0, 

Our  objective  now  is  to  find  the  possible  values  of  e   (i.e., 

the  values  of  e  at  which  buckling  can  start) • 

0  o 

Let  (  )  be  the  value  of  (  )  at  5  =  0,  and  let  (  )  be  the  value 

of  d(  )/d6  at  6  =  0. 

Let  a  =  1  -  e  and  b  =  1  +  (2A8  /X+2u)  so  that  x'  =  ax,,  x»  = 
o  '   o'  1     1   2 

ax^,  and  x'  =  bx-,. 
2*3     3 

Taking  the  doi'ivative  with  respect  to  6  at  0  in  the  equilibrium 
equations  and  boundary  conditions  (I4.03)  -  {h»h)   wo  obtain 


(5.1)    q^i,  1,  =  0(i  =  1,2,3) 

•  •  d  «  c  • 

(5o2)    x'  =  -  Gx,  ,  x'  =  -  exp  ,  and  q,.^  cos  9  +  q,p  sin  Q  =  0  for 

x^  =  R  cos  Q  and  x.-,  =  R  sin  Q'* 
^2  =  0  (i  =  1,2,3)  for  x^ 


(5.3)    qi3  =  0  (i  =  1,2,3)  for  x.  =  -1-  h. 


Next  we  observe  that 


_  '5^ik  _  '^^ik 


'ik,k   3x  ^   op^g  Prs,k  ' 


and  hence 


'lik,k   do^6P^3^ 


'^  ^^ik  '  ''^'^ik  * 

g^Q  Prs,k  "^    ^d9^h=0  Prs,k  =    ^,"^""^6=0  Prs,k 


=    ^°S,^ik'^5=0        ^°^        ^k  =  P,k 


*■,-  ■  '■>■ 


•:-»...'; 
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Also  q^j^  =  O<i^iV^I>^sh=0  ks   =  (^?^li:)5=0  ""^^^^  ^  =  P* 

We  now  modify  (5.1)  -  (5.3)  using  (3.8)  and  the  fact  that  p  is  a 

diagonal  matrix  vrith  positive  diagonal  elements.   Observing  that 

Ap  .  .  -  3A  -  2tx  =  -  2iJ,b  we  obtain  from  (5.1) 

0  =  (Vik^S=0  =  ^^iijj^k  -^  2R;ik,k  -  J-^  ^Pik,k  -  Pki,k)  • 
K  Pii^Pkk 

Prom  (5.2)  when  x^  =  R  cos  Q   and  Xp  =  R  sin  Q  we  have 

0  =  (D^q^^)g^Q  cos  0  +  (0-^32)5=0  ^^^  ^  "^  ^"^^^31  "^  ^Pl3^ 

.  cos  0  +  ^  (ap,2  •*■  ^Ps-^)  si"  ®   » 


and  ■D-,-.  cos  Q  +  p^o  sin  Q  =  0  since  here 


P13  =  ^-  ^^l^x,  =  °   ^"^   P23  =  ^-  '''2^x,  =  ° 


From   (5.3)   when  x.,   =  +  h  we  have 


0  =  (Vl3^5=0  =  'Bk  (^'13  '"  ^P3l^ 


°  =   (^23^5=0  =  Ife  ^'^i^23  ■"  ^^32^      ' 


and 


0  =    {'D^^^^)q=q  =  ?^(Pi3_   +  P22^    "^    ^'^  "^  2[i)p22      . 


Listing   these   results  vie  have 


K  -  r"  ■■ 


{S.h)  Ap.^.^3_   +  ^^"^Vx^lz  - 


2|ab 


p.  .+Pi  , 
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(Pik,k  -  Pki,k^    =  0   (i  =  1,2,3). 


iS.B) 


z 


'    =  -   ex^    ,     x'    =  -   ex^    ,      p,-,    cos    Q  +  p,,^    sin  Q  =  0  for 


X,    =  R   cos   Q  and   Xp   =  R   sin   Q, 


(5.6)  0  =  apj_^  +  ^]?2i  =  aP23  ■*■  ^P^2  ~  ''^^^^ll  "*"  ^22^    +    (A  +  2|a)p 


33 


for  x_  =  +  ho 


To  simplify  the  problem  vie   assume  at  this  point  that  our  solu- 
tion has  radial  syroraetry.   That  is,  we  assume  that  x'  =  f(r,x-)  cos  Q, 

:,  =  r  cos  Q  and  Xp 


x'  =  f(r,x,)  sin  Q,  and  x'  =  g(r,x,)  for  x-,  =  r  cos  Q  and  x^  =  r  sin  Q 


and  for  some  functions  f  and  g, 

Oq  ••  •-  •• 

Then  f  =  ar,   g  =  bx-,   x'    =  f  cos   Q,  x'    =  f   sin  Q,  and   x^,   =  g, 

s>        1  2  :> 

Also 


P  = 


f     cos'^Q  +  =■  f   sin  0 


(f     -  -  f)sin  Q  cos   0 


,    g     cos    Q 
\  ^'r 


(f     -  i  f)sin  Q  cos    0 


f^   sin  e   +  •-  f   cos'^Q 


gp   sin  0 


f     cos    Q 

f     sin  0 
X3 


S,. 


Let  A   =  2(J-a/a+b   and  B  =  A  +    (2iab/a+b)  . 

•  « 

Using  the  fact  that  now  pp-,  =  p,p  vie  obtain  (5.7)  from  (5.i|-)» 


(5.7)   ^ 


(A  +  2m.)(P2i^i  +  P22,2)  "•  ^^^23, 3  "^  ^^32, 3 


,)  +  Ap-j_3^3  -f  Bp33_^3  =  0 

=  0 


(A  4-  2|I)(P3_-L^3_  +  p^2,2 

(A^+  2m.)(P2i^1  +  P22,2'     ^__     ,,^ 
/^P31,l  ■"  ^32,2^  ■"  (^  ■'  2tx)P33^3  +  B(p^^^3  +  P22^3) 


=  0 


Making  full  use  of  radial  symmetry  in  (5^7) j  iB*^)^    and  (5.6)  we 


obtain 


■J  an 


,■.■<:)     0-S^ 


.-  t-       • !  \ 


11 

1  ;.       1 


(5.8) 


]  A(g^^   +  I  S^)    -^    (A   -^  2ix)g  +  B(r^        -.-  i  f      )    =   0 


(5*9)  f (R,x^)    =  -   sR,   g^(R,X2)    =  0 

(5.10)        af^     +  bgj,   =  0    ,    ^(^r   +  p  f )    +    (A  +  2ii)g^     =0     for     ^^  =  +  h. 

•  00 

¥e  look  for  a   series   solution  of    (5.8)    having   the  form  f  =  2;; • 

.        _ci  n-0 

a    (x-)y   (^)    and   g  =  ^ p    (^^-JS    (r)  «     V/hen  we   substitute   these   series 

n  :5  n  ^^_q     n  j  n 

into  (5.8)  we  would  like  to  be  able  to  express  the  contribution  made 
by  the  nth  terras  as  a  product  of  a  function  of  x»  and  a  function  of  r. 
For  this  to  be  possible  it  is  sufficient  that  y'l  +(1/1*) Yl  -  (l/r  )y  » 
Y  ,  and  5'  are  loropertional  and  also  that  S"^  +  (l/'r)Y'>  S  ,  and  y'  + 
(l/r)Yj-  are  proportional.   These  proportionality  conditions  are  satis- 
fied if  Y     =   J^(k  r)  and  5  =  J  (k  r)  vjhere  k  is  an  arbitrary  constant 
and  J  is  the  Bessel  function  of  the  first  kind  of  order  n.   Since  a 
wide  class  of  functions  can  be  expanded  in  series  of  the  functions 
J^(k  r)  or  J  (k  r),  and  since  we  can  easily  satisfy  (5.9)  terra  by  term 
for  such  series,  xve  look  for  solutions  of  the  form 

CO 


f  =  .  er  +2_.a^(^3)Ji(V^ 
n=l     -^ 


00 


and 


S  =  (3,(X3)  +^P,,(x3)J^(k^r) 

vjhere  k  is  chosen  so  that  k  R  is  the  nth  positive  zero  J-,  . 

Using  the  identity  Jq  =  -  J-i  ^^'^    the  differential  equation  for 

J     we  obtain   (d/dr)J   (k  r)    =  -  k  Jt (k  r)    and 
o  '  on  nln 


et-ip.f' 
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+  k  J  (k  r)  •  Viflth  these  results  we  obtain 
'n  o^  n 


and 


Thus 


and 


f^  =  -  e  -.  <^:  a^[.  i  J,(k^r)  -f  ^/^(k^r)]   , 
n=l 


^^   n=l  "  V 


•      co_ 

b;  =  -  3-  k  (3  J,  (k  r )   9 

n^  ^    1^  n  '   ' 


CO 

n=l 


""  '^n^n'?  ^l(  V)  -  Vo'V"   • 


f  +  i  f  =  -  2e  +  '>-  k  G  J  (k  r)   , 

r   r  "^^^^  n  n  o  n    ' 


rr   r  r   ^2       j^.,-  n  n  1'  n  '   ' 


Using  these  results  in  (5.3)  -  (5.10)  x^e  obtain: 

r-(A  +  2ii)k2a^  +  Aa';^  -  Bk^p^  =  0   for   n  =  1,2, 


•  9  • 


(5.11)  I    --^^k^Pn  +  ('^  +  2i.t)P'j!^  +  Bk^^a^!^  =  0   for   n  =  1,2,.,. 
1  p"  =  0 
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aa»  -  bk  p  =  0  (n  =  1,2, .  ..) 


Tor 


(5.12)  <•  Xk^a^   +  (A  +  2tx)(3^I^  =  0  (n  =  1,2,,..) 


x_  =  +  h 
for   x^  =  +  h 


3   - 


-2  As  +  (A  +  a[j,)P'  =  0   for- 


X-  =  +  h 


From  this  we  easily  obtain  (3  =  (2Asx^/A+2ix)  +  (3   for  some  con- 
stant 6   .  Clearly  the  value  of  S   is  determined  by   the  displacement 
■^oo  00  " 

of  the  ^^^hole  plate  parallel  to  the  x_-axis,  and  i-jithout  loss  of  gene- 
rality x^re  may  choose  p   =  0, 

Eliminating  p  from  the  first  two  equations  of  (5»11)  we  obtain 

a^^^  -  2k^a"  +  k^a  =  0  so  that 
n      n  n    n  n 


^n  "  ^'^1^3  "*"  *^2^  ^^^*  ^'n^3  "*"  ^°3^3  "^  ^h)    '^°°^'^   ^'n^'3 


It   then  follows   easily  that 


^n  =  t-   (°3-^3  -"  %^   •*■  ilJ^  "l^"-'"'^  V3  "■  f-    (°1^^3  -^  ^2^   ^  fsf  °3^ 

•    cosh  k  x_      , 


The   first   tx-jo   equations   of    (5.12)    noxj  become 


[(a   +  b)k^(c3X2   +  c^)    +  SB-_bJ_2A+B)_  ^^^^^^^  ^,^^^ 

+   [(a  +  b)k^(c3_X2  +  c^)    +  Siizblf^i^i  c„]cosh  k^x^  =  0 


(5.13)  { 


for  x^  =  +  h-. 


[-2iik^(c^x,   +  c^)    +  ^^{^^)    c^Jsinh  k^x 


3 


+   [-2ixk^(c-x      +  c,  )    +  ^MA±2.l  c^jcosh  k.^x^   =  0  for  x. 


=  +  h. 


V 


Ik 

Me   note  that  in  (5»13)  ^-'e  have  four  linear  homogeneous  equations 

in  the  unknowns  c,  ,c,.,  ,c^,Ci  ,  VJe  make  the  assuraption  that  a  5^  0  and 

S  ^0  for  soiTie  n  >  1  so  that  our  solution  differs  from  the  trivial 
n  — 

solution  to  the  first  order  in  5,   Then  for  some  n  we  have  c,,Cp,c,, 
and  Cj  are  not  all  zero,  and  from  (5»13)  we  obtain  a  determinant  which 
must  be  zero.   Simplifying  that  determinant  equation  we  obtain 

sinh  2k^h        'M2  +  e^)+2[x 

2kir~  "  -  A(2-5sQ)+2r-Tr::2e^   • 

Since  we  ^rant  our  result  to  make  sense  for  plates  v.'hen  e  is 

o 

small,  we  choose  the  plus  sign.   Solving  for  e     we  have 

sinh  2k  h 


2(A+,.)[.— ^-^-  ..  1] 

n ^ 

sinh  2k  X 


(5.11;)   e^        ■  sTnE^iriT 


n 

Let  (T  be  the  Poisson  ratio,  let  D  be  the  depth  or  thickness  of 
the  plate,  and  let  j.  be  the  nth  positive  zero  of  J-,.  Then  (f  = 
X/2{X+[i) ,   D  =  2h,  and  k  =  j'  /Rp  '^^e  can  novj  write  (5.11+)  as 

-r-rr    SipJl   -vr-   -    1 

(5.15)       e^  =         "         " TD     ' 

0+(2+ff-)    Y^  sinJa  -^ 
^n 

Since  (5.11|)  can  be  valid  for  at  most  one  value  of  n  for  each 
value  of  e  ,  our  series  representations  of  f  and  g  become  the  finite 
series 


(5.16)   f  =  -  er  +  a  Jt  (k  r)   ,  g  =  (3  +  p  J  (k  r) 
^    /  n  1  n    '   ^   '  o   "^n  o  n 


where  n  may  be  any  positive  integer,  and  e   is  related  to  n  by  (5. ill-). 
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Using  (5.II4-)  and  (5.13)  and  solvin{^  for  c^jCgjC,,  and  c,  we 


obtain  c- 


0,  C2  = 


wc ,  c^  =  c,  and  Ci  =  Oj,  where  c  is  a  constant 


(not  determined  by  (5<.13)),  and  w  =  h  tanh  l-z^li   +  (aB-b(2A+B))/Bk^(a+b) )  , 
Listing  our  results  for  the  nth  raode  of  buckling  vje  have 


n 


c(-G'slnh  k  x.    +  z_   cosh  k  ^o) 


(5.17)  V  y-[ 


^^'2 


x»    = 


L 


c[-x,   sinh  k  X-,   +    (%Tr—  +Oj)cosh  k-  x-] 
J>  n  _5  n 

ax^  +   [-8X2   +  a.^  -^  J^(k^r)j5  +  0(5^) 
ax2   +   [-ex^   +  a^^  -^  J^(k^r)]5   +  0(6^) 

o 


Now  let  P  be  the  n-th  critical  buckling  pressure.  Clearly 

^  1   «   o    _     1   "   -^ 


p  =  -  11,  so  that  from  (3  =6)  we  have  P^  = -.^   ^-f -jPi -i 

~   "t  ^11*     ^^°^'^   (3.^)   we  have  q^-j_  =   -  2\ib  +  Zas.  so   that 


2:  ^11-^11 


a'-b 


b-a 


2LL(3A+2ii)e, 


Es. 


^n       2(a    .,^         A+2ti+TI^2M--2?:e;re^       l-'^H3(f'l)  ^^'Z^ 


where  E  is   Young's  raodulus   of  elasticity  and  CT"  is    the   Poisson  ratio 
Vje   list   this   as 


(5.18)       P„  = 


^o^ 


"       l-C'+(30'-l)sQ-2o-e2 


where  e   is  given  by  (5»l!4-)  and  (5.15)  . 

It  is  interesting  to  notice  that  e  and  P  both  approach  finite 
limits  (not  dependent  on  n)  as  ^  — >  00.  Hence  for  a  certain  finite 
compression  of  the  cylinder  \-7e   expect  all  modes  of  buckling  to  be 
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solutions   of  the  mathematical  nodel  recardlssn   of   the    thicl-aiess   of   the 
plats, 

6,      Coinparison  x-Jith  Linear  Theor7  Results 
For  thin  plates    e      is   arbitrarily  sraall  and   from   (5.l£)   we  have 

e^E  _co_     l-3J^2cre     ,^  ,, 

^n  -i-.rf^  ^      1-r     ^  'o     • 

Retaining  second  order  terms  in  D/R  in  (5.l5)  we  obtain  e  = 
( 1/12 ( l+C) ) ( j  D/R)  .   Thus  retaining  second  order  terms  in  D/R  in  the 
formula  for  P^  we  have  P^  =  (E/l2(l-o'") )  (  j^D/R)^.   If  vje  let  P  be  the 
flexural  rigidity  of  the  plate,  then  P  =  Eb3/12(1-G^)  and  P^  = 
Pj  /DR  .  Hence  the  nth  critical  buckling  pressure  per  unit  circum- 
ference, P  D,  is  given  by 

(6,1)    P  D  =  -^^- 
n    ^^ 

which  is  the  result  obtained  using  the  linear  theory  for  a  plate  with 
clamped  edges  [3]« 

In  graphs  I  and  II  ue   make  a  graphical  comparison  between  the 
results  of  the  linear  and  nonlinear  theories.   Let  e  ,  be  the  value 
of  e  for  n  =  1,  Graphs  I  and  II  show  how  e  ,  and  P-,/S  respectively 
vary  withD /R  for  the  linear  and  nonlinear  theories  irhen  d' =   0,3. 

If  vje  consider  the  value  of  s  ^  given  by  the  nonlinear  theory  as 
the  correct  value,  and  if  ^^re  think  of  the  value  given  by  the  linear 
theory  as  an  approximate  value,  bhen  for  G^=   Oe3  the  per  cent  error 
due  to  the  linear  theory  is  1,J+  per  cent  for  D/R  =  0,1  and  5.6  per 
cent  for  D/R  =  0,2 «   Thus  the  per  cent  error  Is  small  for  D/R  <  0,1, 
and  it  approaches  zero  as  D/R  — >  0, 
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